Mean-field analysis of the stability of a K-Rb Fermi-Bose mixture 
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We compare the experimental stability diagram of a Fermi-Bose mixture of 40 K and 87 Rb atoms 
with attractive interaction to the predictions of a mean-field theoretical model. We discuss how 
this comparison can be used to give a better estimate of the interspecies scattering length, which is 
currently known from collisional measurements with larger uncertainty. 
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The production of degenerate Fermi-Bose mixtures of 
atoms El [I has opened new directions in the field 
of ultracold atomic gases, with relevant implications for 
the achievement of fermionic superfluidity. Collisional 
interactions between bosons and fermions have a twofold 
role: they allow for an efficient sympathetic cooling of 
the fermions into quantum degeneracy, and they affect 
the properties of the degenerate mixture. In case of a 
large boson-fermion attraction, the most important pre- 
dicted effects are mean-field instabilities PJJE and 
boson- induced interactions between fermions [8j , besides 
a general modification of the properties of the two indi- 
vidual components [iJElEI- We have recently observed 
the collapse of the Fermi gas in a 40 K- 87 Rb mixture [l2| . 
which is indeed characterized by a large attractive inter- 
action between the components. 

In our study of the collapse we found that the order of 
magnitude for the critical atom numbers is in agreement 
with the predictions of a mean-field model for spherical 
geometry 

In this work we investigate the stability diagram of 
the K-Rb system, and we make a quantitative compari- 
son of the experimental findings with mean-field theory 
taking into account the actual trapping potential for the 
atoms. We find that the mean-field model is able to re- 
produce the critical atom numbers for collapse using a 
value for the boson-fermion scattering length a BF that is 
in good agreement with the one we determined through 
the study of cold collisions [Tjj. Moreover, we show how 
the comparison of theory and experiment on the stability 
diagram can be used to determine the scattering length 
with a much lower uncertainty than that determined from 
collisional measurements. 

The Fermi-Bose mixtures of 40 K (fermions) and 87 Rb 
(bosons) atoms are produced by sympathetic cooling in a 
magnetic potential jj| . Both species are trapped in their 
low-field-seeking states with largest magnetic moment, 



which experience potentials of the kind 

Vb(x) = X -ra B [uj 2 b ^{x 2 + y 2 ) + u? Bz z 2 \ (1) 

Vf(x) = ^m F [u; 2 F± {x 2 + {y-y a ) 2 )+uj Fz (z- z ) 2 ] . 

The trap frequencies are co F j_ = 2tt x 317 Hz, uj Fz = 
2i r x 24 H z for 40 K, while those of 87 Rb are a factor 
\J m B jm F ~ 1.47 smaller. The system is characterized 
by a gravitational sag between the two clouds, which af- 
fects both the horizontal and vertical directions, due to a 
small misalignment of the magnetic trap with respect to 
the direction of gravity 01 ■ The values of the displace- 
ment between the two potential centers are yo = 3.6 /zm 
and zq = 20 /im. 

The largest samples that we are able to cool to quan- 
tum degeneracy amount to about 2xl0 5 bosons coexist- 
ing with 5x10 fermions. The atom numbers for both 
species can be adjusted separately by controlling the 
loading and the evaporation procedures. The coldest 
Fermi gas observed have a temperature of about 0.2T F , 
and the minimum temperature of the bosons that we are 
able to measure from the non-condensed fraction is lim- 
ited to about 0.6T C . 

In the mean-field model that we use to describe this 
Fermi-Bose system, the ground state is determined by 
the solution of Gross-Pitaevskii (GP) equation for the 
bosons, coupled to the Thomas-Fermi equation for the 
fermions |5j,|6J 
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where n B — \<t> 2 \ and np are the boson and fermion 
densities respectively, and </»(x) is the condensate wave 
function. The boson-boson and boson-fermion interac- 
tions are described by the coupling constants g BB = 
^■nh 2 a BB /m B and g BF = 2nTi 2 a B p /m,R in terms of the 
s-wave scattering length a BB and a BF . Here m B ^ F are 
the atomic masses and mji — m B m F /(m B + m F ) is the 
reduced mass. 
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In the case of the rubidium-potassium mixture, the 
best experimental estimates of the scattering lengths are 
a BB = 98.98 ±0.02 a [3 and a BF = -410 ±80 a [3. 
The difference in the uncertainties for the two quantities 
stems from the different techniques used, which is a com- 
bination of photoassociation and Fcshbach spectroscopies 
in the first case, and measurement of cold collisions in the 
second case. Although the investigation of cold collisions 
is the easiest way to measure the interactions in an ul- 
tracold gas, it is often accompanied by a rather large 
uncertainty. Indeed, what one measures in the experi- 
ment is a collisional rate, which at ultralow temperature 
has a dependence 

T oc No 2 (4) 

on the atom number N and on the scattering length a. 
The absolute atom number, in the case of a mixture of 
two species, cannot be determined with an accuracy typ- 
ically better that 40%, giving rise to a 20% uncertainty 
on the scattering length. 

By solving equations ©-© we can investigate the 
ground state and the stability of the system against col- 
lapse. The equations are solved with an imaginary-time 
evolution, embedded in an iterative scheme, as in [(|. 
The stability is checked by requiring that energies, chem- 
ical potentials and central densities converge to the final 
value, with a precision at least of 10 -6 . In Fig. ^ we show 
a typical ground state configuration, calculated with the 
nominal values a B F =-410 a , Nr = 7.5-10 4 , N F = 2-10 4 . 
As already discussed in Refs. the mutual attrac- 

tion results in an enhancement of the density of both 
species in the volume of overlap. We note that the grav- 
itational sag, that shifts the centers of the two clouds, 
reduces the effect of the mutual attraction with respect 
to the concentric case studied in Ref. . 

The deformed ground state of the system is predicted 
to lead to a series of effects that can be described within 
a mean-field approach, including a modification of the 
frequencies of collective excitations |9j, [llj and of the ex- 
pansion of the two clouds from the trap [HJ- We ac- 
tually observed the latter effect already from the first 
experiments 4] , where we had an evidence of a modified 
expansion of the condensate in presence of the Fermi gas. 

If the atom numbers are increased above some critical 
values, an instability occurs. In our model, accordingly 
to the study reported in Refs. 0, 0] the signature of 
the instability is the failure of the convergence procedure 
during the iterative evolution toward the ground state 
of the system. In particular, the onset of instability is 
characterized by an indefinite growth of central densi- 
ties which triggers the simultaneous collapse of the two 
species. We note that the investigation of the actual dy- 
namics after the system has been driven into the unstable 
region would require a description that goes beyond Eqs. 
iJUl-©, in a similar fashion to what happens during the 
collapse of a single Bose-Einstein condensate with attrac- 
tive interaction 0, 0| . Here we will not discuss these 
aspects, since we are concerned with the determination 
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FIG. 1: Density profiles along the vertical y direction (a) and 
contour levels in the (y — z) plane (b) of the ground-state of 
the Fermi-Bose mixture, for a_BF=-410 ao, Nb = 7.5 ■ 10 4 , 
N F = 2 • 10 4 . 



of the critical values for the onset of instabilities only. 

In the experiment we observe the instability of the 
Fermi-Bose system as the number of condensed bosons in 
increased by evaporation, at a fixed number of fermions 
[l2| . The signature of the instability is a sudden loss of 
most of the fermions from the magnetic trap, which we 
attribute to a largely increased three-body recombination 
in the collapsing samples. Here we are interested just in 
the critical atom numbers, but we note that further ex- 
periments, for example using a Feshbach resonance to 
tune a B F |13| . will be necessary to study the collapse 
dynamics of the system. 

To compare the predictions of the mean-field model 
to the experimental findings on the instability, we have 
built a stability diagram, shown in Fig. [21 Here we re- 
port in the plane N B ~ Np the condensate and fermion 
atom numbers that we were able to measure in the ex- 
periment for stable samples, and compare them with the 
calculated critical line for instability. We note that all 
the data points refer to samples at temperatures of the 
order of 0.2-0.5Tp for the Fermi gas, and with almost no 
detectable thermal fraction for the BEC. In the model, 
the occurrence of instability depends on both N B , N F , 
and also on the value of the scattering length a BF . In 
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FIG. 2: Region of stability of the Fermi-Bose mixture, as a 
function of the number of atoms. The black dots are the 
experimental points; lines are the theoretical prediction for 
the boundary between the stable (left) and collapse (right) 
regions, for three values of the inter-particle scattering length 
(in units of the Bohr radius ao): clbf — —380 ao (dotted line), 
aBF = —395 ao (continuous line), clbf = —410 ao (dashed- 
dotted line). The marked dots are found very close to the 
instability (see text). 

Fig. 0we have therefore plotted a family of critical lines 
by varying clef around the nominal value. Note that the 
position of the critical line depends quite strongly on the 
value of aBF, and it shifts toward larger atom numbers 
for decreasing magnitudes of clef- Although the critical 
line calculated for the nominal value asF=-410 ao is m 
rather good agreement with the experimental observa- 
tion, the one calculated for clef =-395 ao better fits the 
data. In the experiment we observed the collapse of the 
Fermi gas for number pairs close to the two marked data 
points (l9j, which appears to be actually the region of 
achievable number pairs closest to the instability. 

The strong dependence on the criticality on asF that 
we find in the mean-field analysis is in accordance with 
the scaling law for the critical number of condensate 
atoms 

N c B Ht ~\a BF \- a (5) 

with a = 12, discussed in [5j,[3|. Here we have verified 
that this scaling behavior perfectly fits the numerical so- 
lutions of Eqs. JU-JU! in case of a simplified concentric 
spherical configuration, by using an average trapping fre- 
quency as usually found in literature 0, l2f)j | , as shown in 
Fig. |3 In this calculation we have kept the number 
ratio Ns/Np fixed, which corresponds to move along a 
straight line passing from the origin in the diagram of 
Fig. This calculation also shows that the 3D geome- 
try of the experiment results in a renormalization of the 
scaling exponent, a^,D — H, besides an effective increase 
(in modulus) of the critical scattering length, of the or- 
der 2%, due to the presence of gravity. In fact, the effect 
of the gravitational shift is a reduction of the density 



FIG. 3: Behavior of the critical number of atoms Nb as a 
function of the scattering length clef, for Nb = 4.25 Nf (tri- 
angles) and Nb = 8.5 Nf (circles). The numerical solutions 
of the mean-field equations (symbols) are compared with the 
scaling law in Eq. JSJ with exponent a = 12 (continuous 
line) and a = 11 (dashed line). Filled and empty symbols re- 
fer respectively to the full 3D geometry and to the simplified 
spherically symmetric case discussed in the text. 

enhancements due to the attractive interaction, with re- 
spect to the case in which the to species are trapped by 
concentric potentials. 

This analysis suggests the possibility of using a com- 
parison of the experimental and calculated critical atom 
numbers to extract a precise information on the value of 
the scattering length for our K-Rb mixture. To account 
for the uncertainty in the determination of atom num- 
bers in the experiment, in Fig. 0] we have plotted the 
same data points of Fig. |2 expanding or compressing 
both scales by 40%, and compared them with the same 
family of critical lines. 

By taking the critical lines that best fit with the ob- 
servation in these two limiting cases, we can extract the 
following mean-field estimate for gbf 

a% ~ -395 ± 15 a . (6) 

Note how the dependence of as f on the number of atoms 
of Eq. [SJ stronger than that of Eq. results in a 
5-fold smaller uncertainty on the scattering length de- 
termined by the stability analysis than that obtained 
with collisional measurements. We expect this mean- 
field prediction for obf in Eq. © to be rather robust, 
since the 40% experimental uncertainty on the number 
of atoms should be quite larger that the possible discrep- 
ancy between the GP theory and the actual behavior of 
the system, as recently discussed for the case of a single 
85 Rb Bose-Einstein condensate with attractive interac- 
tion [2ll2ll23|. 

The assumption of the same scaling for boson and 
fcrmion atom numbers due to experimental uncertainty 
is not unrealistic. Indeed, the possible sources of sys- 
tematic errors in the atom number calibration are resid- 
ual magnetic fields, non perfect polarization of the probe 
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FIG. 4: Region of stability of the Fermi-Bose mixture, as 
a function of the number of atoms. To keep into account 
the experimental uncertainty on the number of atoms these 
have been increased or reduced by a 40%, respectively in (a) 
and (b). Continuous lines are the theoretical prediction for 
the boundary between the stable (left) and collapse (right) 
regions, for three values of the inter-particle scattering length 
(in units of the Bohr radius ao): cibf — —380 ao (dotted line), 
cibf = —395 ao (continuous line), cibf = —410 ao (dashed- 
dotted line). 



beam, or improper calibration of the imaging magnifica- 
tion. When using a simultaneous imaging scheme |J] , all 
these factors affect in the same way the measured atom 
numbers for bosons and fermions. 

We note that recently it has been proposed a beyond 
mean-field approach, including the second order correc- 
tion in the scattering length cibf, that should provide 
important modifications to the mean- field results |20| . 
Nonetheless we have verified that the inclusion of this 
term stabilizes the system in the range of atom numbers 
considered here for any value of the scattering length, 
thus forbidding the collapse, in contrast with the ex- 
perimental findings. Higher order terms might therefore 
be important to provide a correct beyond mean-field de- 
scription, although at this stage it is not clear whether 
this could affect the critical number of atoms for collapse. 

In conclusion, our mean-field analysis of the stability 
of the K-Rb Fermi-Bose mixture shows the possibility 
of an independent determination of the relevant scatter- 
ing length of the system. The strong dependence of the 
critical atom numbers on the scattering length allows to 
largely reduce the effect of the experimental uncertainty. 
The value we determine for clbf is in very good accor- 
dance with that found from collisional measurements. 
We note that an independent, more precise determina- 
tion of clbf, for example using Feshbach spectroscopy 
(l3| , would be useful to test the validity of the mean-field 
theory for strongly interacting Fermi-Bose mixtures, and 
to asses the importance of beyond mean-field effects. 
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